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ABSTRACT 

Current  meter  data  were  collected  from  moorings  on  the  800  and  1800  m  isobaths 
on  the  continental  slope  off  Point  Sur,  California  at  100,  350,  and  500  m  depth  and  at 
1000  m  depth  on  the  1800  m  isobath.  Complex  demodulation  and  spectral  techniques 
were  used  to  study  the  distribution  of  tidal  energy  at  the  diurnal  and  semidiurnal  fre- 
quencies, and  to  study  the  size,  rotation  and  orientation  of  the  tidal  current  ellipses.  At 
mooring  P2  (800  m)  the  observed  diurnal  energy  increased  with  depth  for  each  of  the  17 
months  of  data  analyzed,  and  the  semidiurnal  energy  likewise  increased  with  depth  for 
14  of  the  17  months.  The  observations  were  studied  using  the  theory  of  bottom  trapped 
waves  in  a  rotating  stratified  fluid  (Rhines,  1970).  Bottom  trapping  is  expected  at  the 
diurnal  and  semidiurnal  frequencies  whenever  the  critical  trapping  frequency 
coe  =  /Vsinasin</>  exceeds  1/24  and  1/12  hours-1  respectively,  where  N  =  the  Brunt- Vaisala 
frequency,  a  the  bottom  slope,  and  </>  is  a  measure  of  the  angle  the  wavenumber  vector 
makes  with  the  gradient  of  topography.  Theoretical  energy  decay  curves  matched  the 
observed  energy  levels  at  P2  very  well,  indicating  that  the  enhancement  of  energy  with 
depth  was  likely  due  to  the  presence  of  bottom  trapped  waves.  No  bottom  trapping  was 
observed  at  P3  (1800  m)  because  the  deepest  current  meters  were  positioned  too  far 
above  the  sea  floor  to  detect  bottom  trapped  energy,  and  because  the  density  gradients 
near  the  sea  floor  were  weaker  than  at  site  P2  and  kept  we  below  the  tidal  frequencies. 
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I.     INTRODUCTION 

A.     BACKGROUND 

The  Naval  Postgraduate  School  maintained  an  array  of  subsurface  current  meter 
moorings  ofT  Point  Sur,  California  from  1989  to  1990.  The  moorings  were  positioned 
on  the  continental  slope  and  collected  data  used  in  the  study  of  the  California  Current 
System  (Figure  1).  The  initial  spectral  analysis  of  data  from  mooring  P2  indicated  that 
the  kinetic  energy  at  diurnal  tidal  frequencies  increased  with  depth  (Figure  2).  This  was 
an  unusual  phenomena  because  vertical  energy  structures  at  tidal  frequencies  are  usually 
associated  with  the  baroclinic  modes  of  free  internal  waves.  At  the  latitude  of  Point  Sur 
(36°  20'  N)  the  inertial  frequency  is  approximately  1.1  cpd  and  free  internal  waves  can- 
not exist  at  diurnal  (subinertial)  frequencies.  Baroclinic  waves  can  exist  at  subinertial 
frequencies  only  if  the  waves  are  trapped  along  a  wall  or  slope. 

The  tides  off  most  of  the  California  coast  are  a  mixture  of  diurnal  and  semidiurnal 
components.  The  principal  lunar  constituent,  M2,  is  the  strongest  tidal  element  followed 
by  the  Ku  Ou  and  S2  constituents.  Barotropic  diurnal  and  semidiurnal  tides  both  propa- 
gate northward  along  the  coast.  The  semidiurnal  tide  is  a  combination  of  Kelvin, 
Poincare  and  forced  waves.  The  barotropic  semidiurnal  tidal  current  ellipses,  represent- 
ing horizontal  current  amplitudes  and  directions,  are  generally  oriented  parallel  to  the 
coastline  (Munk  el  a/.,  1970).  The  barotropic  diurnal  tide  is  a  combined  Kelvin  and  shelf 
wave  with  current  ellipses  which  can  be  either  parallel  or  perpendicular  to  topography 
(Noble  ei  al,  1987).  Several  studies  of  California  tides  have  found  the  primary  direction 
of  rotation  for  diurnal  and  semidiurnal  current  ellipses  to  be  counterclockwise. 

Baroclinic  tides  are  generated  when  semidiurnal  surface  tides  force  isopycnal  sur- 
faces over  topographic  features,  such  as  the  continental  shelf  break  (IVunsh,  1975). 
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Figure   1.       Point  Sur  Moorings  P2  and  P3:     The  location  of  Point  Sur  on  the 
California  coast  is  shown  in  the  inset. 


While  the  barotropic  semidiurnal  tide  is  usually  aligned  with  topography,  internal  tides 
often  become  the  dominant  cross-slope  current  component  along  the  California  coast 
(Noble  el  ai,  1987). 

Variations  in  the  vertical  distribution  of  tidal  energy  have  been  observed  at 
semidiurnal  frequencies.  During  the  Coastal  Ocean  Dynamics  Experiment  (CODE),  the 
baroclinic  semidiurnal  internal  tide  was  larger  than  the  barotropic  tide.  Changes  in 
semidiurnal  kinetic  energy  at  different  depths  were  linked  to  surface  winds  and  changing 
stratification  (Roscnfcld,  1990).  The  semidiurnal  internal  tides  observed  in  CODE  had 
some  bottom  energy  intensification  which  was  attributed  to  the  horizontal  velocity 
structure  of  the  first  dvnamic  mode. 
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Figure  2.       Spectral  Analysis  of  Energy  at  Mooring  P2  During  July  1989:     The 

curves  represent  energy  at  100  m,  350  m  and  500m. 


There  are  few  published  studies  which  discuss  similar  bottom  intensification  of 
subinertial  currents.     Bottom  trapped  topographic  Rossby  waves  were  observed  by 


Thompson  and  Luyten  (1976)  and  were  discussed  using  bottom  trapped  wave  theory 
developed  by  Rhines  (1970).  These  were  low  frequency  waves,  however  similar  theory 
was  applied  to  bottom  trapped  waves  at  tidal  frequencies  and  used  to  study  the  tidal 
currents  and  vertical  energy  distribution  of!  Point  Sur. 

B.     BOTTOM  TRAPPED  WAVES  IN  A  ROTATING  STRATIFIED  FLUID 

This  discussion  follows  the  development  presented  by  Rhines  (1970)  for  bottom 
trapped  waves  in  a  uniformly  stratified,  rotating  fluid  bounded  by  a  single  plane  wall 
inclined  everywhere  at  angle  a  to  the  horizontal.  A  right  handed  orthonormal  coordi- 
nate system  is  oriented  on  the  slope  such  that  the  jr,  and  x2  axis  are  rotated  at  an  angle 
4>  from  the  alongslope  and  upslope  directions  respectively  and  x3  is  perpendicular  to  the 
bottom  slope  (Figure  3).  The  corresponding  velocities  along  xu  x2  and  x3  are  designated 
as  w,  ,  w2  and  u3  . 

The  governing  equations  are  the  linearized,  non-diffusive  Boussinesq  equations  used 
in  the  form: 


^  +  2Qu  =  -^-,p'f  (1) 


dp'        A'2 


di         S 


uj  (2) 


where  p  =  pressure,  Q,  =  angular  rotation  of  the  earth,  g  =  gravitational  acceleration, 

j     =    a  vertical  unit   vector,   2Q.  =   f\    =    f(  sin  a  sin  (/>,     sin  a  cos  0,     cos  a)   ,  /= 

—g   dpx 
Coriolis  parameter,  the  Brunt-Vaisala  frequency  squared  =  N2  =  — —  ,  and  a  per- 
turbation density  is  used  such  that  p  =  p0  +  p(z)  +  p'{x,  y,  z,  t)  . 

Inviscid  boundary  conditions  are  chosen  so  that  w3  =  0  at  jc3»0.    If  u,»0  every- 
where in  the  water  column  the  motion  is  nondivergent  in  planes  perpendicular  to  the  x3 

dux        du2 
axis  and  the  continuity  equation  — —  +  — —  =  0  applies.    The  wave  propagates  in  the 

CX\  CX-) 


Figure  3.      Coordinates  for  Bottom  Trapped  Waves  on  a  Sloping  Bottom:     The  ele- 
vation from  the  horizontal  is  given  by  a. 


du. 


du, 


x2  direction  with  — —  =  0  .    To  balance  the  continuity  equation  -—-  must  also  equal 
oxx  ox2 

a  2 

zero.    When  u2  j*  0  and  — —  =  0,  free  internal  waves  can  exist.    If  u2  =  0  and  — —  #  0 

ox2  ox2 

waves  are  trapped  against  the  wall  (bottom)  in  a  plane  parallel  to  the  boundary  surface 
with  particle  motion  limited  to  the  jc,  direction. 

Substituting  ult  u2  =  u3  =  0,  and  the  appropriate  form  of  2Q,  into  Equations  (1)  and 
(2)  leads  to  a  set  of  three  equations  where  for  the  x}  direction,  Equation  (1)  becomes 


ou] 
-T-  =  -gp'(  sin  a  sin  <f>)  (3) 

01 


In  the  x2  direction,  Equation  (1)  becomes 


1    /   dp         dp         dp 
/u,(coS«)  =  -  —  /—  +  —  +  —  I 


1      °p 
/w,(cosa)  =  -  — gp'{  sin  a  cos  0)  (4) 


For  the  ar3  direction,  the  result  is 


-fU]{  sin  a  cos  4>)  =  -~pj—-  SP'(  cos  a)  (5) 


The  rate  of  change  of  the  density  perturbation,  equation  (2),  similarly  reduces  to 


dp'        v2 

— - —  =  •±rr~  w,(  sin  a  sin  <b).  (6) 

ct         S 


The  combination  of  equations  (3)  and  (6)  results  in  an  expression  for  the  frequency  a 
where: 


dp'  \(  du, 

»  =  (iV2t/,p')(  sin  a  sin  4>)2 


dt      \    dt 


-.2 

-~  (p'Mj)  =  N2U]p'(  sin  a  sin  4>f 
dt 


a2  =  N2(  sin  a  sin  </>)2 


a  =  ±  N  sin  a  sin  <fi  (7) 


Substituting  (7)  into  (6),  p'  can  be  determined  as 


p'— J-"i  <8) 


Solutions  to  the  wave  equations  are  found  using  the  form 

p=pQe(-K  +  *>)*  +  W2*2-'")  (9) 

where  k  and  A3  are  real,  k  is  defined  as  a  trapping  coefficient  with  units  of  wavenumber. 
The  alongslope  or  horizontal  (  k2  )  and  vertical  (  A:3  )  wavenumbers  in  this  solution  lead 
to  spatial  dependence  in  the  x2  and  jc3  directions.  Using  equations  (8)  and  (9)  in  (4)  and 
(5),  the  spatial  dependence  of  the  wave  is: 

k  5[1  -  sin2a  sin2(ft] 

— —  sgn(  sin  a  sin  <p)  (10) 


c2        cos  a  +  S   sin  a  cos  0 

k,        1      (1  —  S  )  sin  2a  cos  6 

— -  =  —  — —  (11) 

K2       L    cos  &.  +  S   sin  a  cos  4> 

N 
where  S  =  — .    When  the  trapping  coefficient  k  is  positive  the  wave  motion  is  expo- 
nential and  trapped  near  the  bottom.   From  equation  (9),  wave  motion  decreases  due  to 
the  effect  of  trapping  by 

g"**3  (12) 

The  magnitude  of  trapping  is  dependent  upon  a  and  </>.  When  particle  motion  is  di- 
rectly up  the  slope  {<fi  =  90°)  trapping  is  maximized.  When  particle  motion  is  along  the 
isobaths  (</>  =  0°),  no  trapping  occurs. 

Rhines  discusses  four  different  geometries  where  trapping  occurs:  (1)  a  ~  1  for  steep 
slopes,  (2)  a  ~  S_l  for  more  gradual  slopes  such  as  those  on  the  continental  slope,  (3) 


a<^5"'  for  very  gradual  slopes  such  as  those  on  the  continental  rise,  and  (4) 
—  —  a  <  S  cos  (f)  for  near  vertical  slopes.  Consider  the  case  appropriate  for  the  conti- 
nental slope  off  the  California  coast  where  a  ~  S"1  and  sin  c/>  ~  O(l).  Evaluating 
equation  (7)  using  these  relationships  leads  to  the  following: 


a    =     ±N  sin  a  sin  0    ~    A"  sin(S~  )  sin  <$>    ~    N-r? 

A 


/  d3) 


— —  can  be  expressed  as: 


K 


S[_  1  —  sin'S     sin  (/>] 


c2  co^S     +  5'  sin  S     cos  4> 


-    S  (14) 


The  ratio  of-2-  reduces  to 


_*3_  5[1  -sin2(S~')sin  </>] 

ki  cos:(^_1)  +  S2  sin2(5"!)  cos20 


=    5  (15) 


The   relationships   now  established   with      a  -*  Sl      arc      a  ~f,       — S,      and 

"7"^  S.  Parallel  to  the  slope  Coriolis  is  balanced  by  buoyancy  as  shown  in  equation 
(15).  Coriolis  traps  waves  against  the  wall  (bottom)  but  the  frequency  of  these  trapped 
waves,  determined  by  equation  (7),  is  not  dependent  on  /  Bottom  trapped  internal 
waves  are  similar  to  internal  Kelvin  waves  because  they  are  not  bound  by  a  lower  fre- 
quency limit  of  /and  trapped  waves  can  exist  and  propagate  at  subinertial  frequencies. 
The  high  frequency  limit  of  bottom  trapped  waves  is  a  and  is  a  function  of  a  and  N. 
Maximum  trapping  will  always  occur  when  </>  =  90°  and  the  particle  motion  is  directly 
up  the  slope. 
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The  exponential  decay  of  trapped  waves  away  from  the  bottom  is  determined  by 
trapping  coefficients.  To  determine  k,  bottom  trapped  waves  were  assumed  to  have  a 
dispersion  relationship  similar  to  baroclinic  Kelvin  waves  (Rhines,  1970)  where 


v2,  2 

n=  1,2,  3  ... 


*Mit-) 


with  a)  the  frequency  of  a  wave,  k  a  horizontal  wave  number,  n  the  mode  number  and 
H  the  water  depth.  The  use  of  this  dispersion  relationship  represents  an  approximation 
because  baroclinic  Kelvin  waves  assume  a  flat  ocean  bottom.  The  horizontal  wave 
number  can  then  determined  using 


,2         2    V   H 

k   =  co 


(it)2 


"    /V2  2s 

(A    -oyn) 


When  <f>  =  90°  waves  are  propagating  along  the  slope  in  the  x2  direction  and  particle 
motion  is  up  and  down  the  slope  in  the  jc,  direction.  Trapping  coefficients  are  defined 
as  k  =  k2S  from  equation  (14)  and  can  be  calculated  using 


7(  nn  \; 

wAir) 


/  V    (N2  -  col) 


Waves  are  trapped  closer  to  the  sea  floor  and  decay  more  rapidly  away  from  the 
bottom  when  k  increases.  Trapping  coefficients  increase  when  (1)  the  buoyancy  fre- 
quency AT  increases,  (2)  depth  H  decreases,  (3)  higher  frequencies  o>  are  considered  (in 
the  tidal  range  of  1.0  to  2.0  cpd),  (4)  f  decreases  (lower  latitudes)  and  (5)  for  higher 
modes  n.  Trapping  is  weakest  and  waves  decay  more  gradually  away  from  the  bottom 
in  the  first  baroclinic  mode.   No  trapping  occurs  for  the  barotropic  case  when  n  =  0. 


II.     DATA  COLLECTION 

A.     CURRENT  METER  DATA 

Current  data  were  collected  at  two  sites  ofT  Point  Sur,  California  (Tabic  1)  located 
over  the  continental  slope  (a  =  .06). 

Table   1.     POINT  SUR  MOORING  LOCATIONS 


Mooring 

Position 

Depth 

P2 

36-20.0  N    122-10.2  W 

800  m 

P3 

36-20.0  N    122-27.6  W 

1 800  m 

Mooring  P2,  on  the  800  m  isobath,  was  equipped  with  three  Aanderaa  vector  averaging 
current  meters  (RCM8)  at  depths  of  100  m,  350  m  and  500  m.  Mooring  P3  was  situated 
52  km  offshore  in  1800  m  with  four  Aanderaa  current  meters  at  depths  of  100  m,  350 
m.  500  m  and  1000  m  (figure  4). 

Vector  averaged  current  velocities,  current  directions,  temperature,  conductivity  and 
pressure  measurements  were  recorded  every  30  minutes  in  a  solid  state  Data  Storage 
Unit  (DSL).  The  DSU  had  a  storage  capacity  of  65530  10-bit  words  and  could  record 
seven  months  of  data  with  the  sampling  interval  used.  Sampling  times  were  regulated 
by  an  internal  quartz  clock  with  an  accuracy  of  ±2  seconds  per  day.  Several  deploy- 
ments were  necessitated  by  the  seven  month  data  storage  limit  (Figure  5).  P2  was  re- 
covered and  redeployed  in  May,  August  and  December  1989  and  May  1990,  The  three 
current  meters  used  on  P2  provided  nearly  continuous  data  from  May  1989  through 
October  1990.  Mooring  P3  was  deployed  in  December  1989  and  again  in  May  1990. 
Good  data  was  obtained  at  four  levels  from  December  1989  to  April  1990.  After  April, 
the  data  from  the  current  meters  at  500  m  and  1000  m  was  limited  by  instrument  failures. 
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Figure  4.      Current  Meter  Configuration  at  Moorings  P2  and  P3 

Velocities  were  measured  using  a  shrouded  paddlewheel  at  the  top  of  the  recording 
unit.  The  paddlewheel  was  magnetically  coupled  to  an  electronic  counter  and  velocities, 
in  units  of  cm  sec1,  were  determined  by  the  number  of  revolutions  the  paddlewheel 
made  during  each  sampling  interval.  The  sensor  had  a  range  of  2  -  250  cm  sec  '  with 
an  accuracy  of  ±  1.0  cm  sec-1  or  +  2  percent  of  the  current  velocity.  The  paddlewheel 
stalled  at  velocities  below  a  threshold  of  2.0  cm  sec  '  (Aanderaa,  1990).  The  cosine  re- 
sponse was  acceptable  within  ±  27°  from  the  vertical.    Both  the  pressure  sensors  and  the 
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Figure  5.      Current  Meter  Deployments  at  Moorings  P2  and  P3:     The  shaded  bars 
indicate  the  intervals  where  good  data  was  obtained  at  each  depth. 

mooring  dynamics  program  in  use  at  NPS  indicated  a  maximum  mooring  tilt  of  11.7° 
which  was  well  within  acceptable  limits. 

The  current  meters  were  oriented  into  the  current  by  a  large  vane  with  an  area  of 
2425  cm2.  The  direction  sensor  was  a  magnetic  compass  with  a  potentiometer  located 
at  the  bottom  of  the  recording  unit.  The  compass  was  oil  damped  and  required  3  to  5 
seconds  to  resolve  changes  of  90°  or  more.  The  resolution  of  the  compass  was  0.35° 
with  an  accuracy  of  ±  5°  for  velocities  5  -  100  cm  sec-1  and  ±  7.5  degrees  for  velocities 
of  2.5  -  5.0  cm  sec-1  and  100  -  200  cm  sec-1.  The  compass  could  function  properly  with 
a  maximum  tilt  angle  of  12°  from  the  vertical.  Prior  to  each  deployment  the  combined 
magnetic  deviation  and  variation  of  each  unit  was  determined  on  a  surveyed  test  bench 
at  the  Naval  Postgraduate  School  and  the  necessary  corrections  were  applied  via 
postprocessing  software. 
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B.     HVDROGRAPHIC  DATA 

While  the  current  meters  were  deployed,  several  oceanographic  research  cruises  were 
conducted  in  the  region  by  the  Naval  Postgraduate  School.  These  voyages  consisted  of 
California  Undercurrent  (CU)  and  Student  Cruises  (ST)  which  collected  CTD  data  at 
intervals  of  approximately  10  km  along  the  Point  Sur  Transect  (POST).  Data  were  col- 
lected near  the  moorings  from  the  surface  to  within  50  -  150  m  of  the  bottom  using  a 
Neil  Brown  Instruments  System  (NBIS)  Mark  1 1  IB  CTD  and  are  described  in  detail  in 
Tisch  (1990).  Using  the  data  from  the  cruises  (Table  2),  the  various  characteristics  of 
the  waters  off  Point  Sur  were  examined.  Although  the  station  locations  and  spacing 
varied  slightly  with  each  cruise,  the  CTD  data  necessary  for  interpretation  of  the  current 
meter  results  were  comparable. 


Table  2.     HYDROGRAPHIC   CRUISES   OCCURRING   DURING 
CURRENT  METER  DEPLOYMENTS 


Cruise  Designation 

Date  of  CTD  Casts 

STMAY89 

5/5/89 

CUJUL89 

7/29/89 

CUSEPS9 

9/26/89 

STNOV89 

11/2/89 

CUJAN90 

1/18/90 

CUMAR90 

3 '2 '90 

STMAY90 

5:2  90 

CUJUN90 

6'27'90 
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III.     METHODS  OF  ANALYSIS 

A.     SMOOTHING  AND  FILTERING 

Upon  recover}',  the  current  meter  DSUs  were  downloaded  into  a  computer.  The  raw 
data  files  were  converted  standard  scientific  units  using  the  sensor  calibration  equations 
and  stored  as  formatted  ASCII  files.  All  data  records  were  plotted  to  facilitate  visual 
inspections.  Velocity  records  were  truncated  where  instrument  failures  or  malfunctions 
were  evident.  Spurious  values,  such  as  an  occasional  outlier,  were  edited  or  deleted 
manually. 

All  data  were  band  pass  filtered  prior  to  further  analysis  to  isolate  the  frequencies 
of  interest.  The  filtering  and  interpolation  methods  used  were  originally  developed  by 
Oregon  State  University  (Denbo  et  al.  1984).  Postprocessed  records,  containing  data 
at  30  minute  intervals,  were  first  low  passed  filtered  (LP)  using  a  Cosine-Lanczos  filter. 
Such  a  filter  minimizes  unwanted  ringing  or  oscillations  caused  by  the  Gibbs  phenomena 
while  still  providing  a  sharp  roll  off  at  the  desired  cutofT  frequency.  A  centered  25  point 
segment  of  the  records  was  used  to  produce  each  LP  filtered  data  point.  The  LP  filter 
had  a  half  power  period  of  2.9  hours  (or  a  frequency  of  8.4  cpd).  The  filter  allowed  any 
signal  with  a  period  exceeding  3.6  hours  to  pass  through  unattenuated.  Signals  with  a 
period  of  less  than  2.0  hours  were  removed  entirely.  After  the  LP  filtering,  records  were 
interpolated  to  transform  30  minute  data  records  to  60  minute  intervals  with  each  data 
point  falling  on  an  even  hour.  The  interpolation  routine  utilized  Lagrangian 
polynomials  to  interpolate  any  unevenly  spaced  data  points  (due  to  missing  individual 
data  points  or  outliers)  to  specific  time  intervals  {Gerald  and  IVheatly,  1989).  After  in- 
terpolation, the  data  was  again  low  passed  filtered  (LLP).  The  second  Cosine-Lanczos 
filter  had  a  half  power  period  of  46.59  hours  (or  a  frequency  of  0.515  cpd)  and  used  a 
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centered  121  point  segment  to  produce  one  filtered  data  point.  This  filter  was  designed 
to  remove  diurnal  and  shorter  period  tidal  energy  and  any  inertial  frequency  energy 
present  in  the  signal.  The  LLP  records  were  subtracted  from  the  LP  records  to  produce 
the  band  passed  form  used  in  the  tidal  analysis.  A  sample  from  100  m  at  P2  (Figure  6) 
shows  the  results  of  the  filtering  procedures. 

B.    TIDAL  COMPONENTS 
1.     Harmonic  Analysis 

To  understand  the  nature  of  the  tidal  currents,  an  analysis  of  the  various  tidal 
components  was  desired.  Traditional  harmonic  tidal  analysis,  or  the  technique  of  linear 
combinations,  requires  a  complete  set  of  observations  over  a  standard  length  period. 
Tidal  heights  or  motions  can  be  expressed  as  a  sum  of  harmonic  terms,  each  determined 
by  a  known,  fixed  frequency.  A  single  estimate  of  the  amplitude  and  phase  are  then 
extracted  through  the  analysis  of  observations.  Observational  periods  are  normally- 
multiples  of  synodic  periods  involving  the  most  important  tidal  constituents.  The  pre- 
ferred record  length  for  a  complete  tidal  analysis  is  369  days.  This  allows  extraction  of 
the  amplitude  and  phase  of  at  least  64  of  the  largest  tidal  constituents  (Dronkers,  1964). 
An  adequate  description  of  the  major  constituents  can  be  accomplished  using  a  much 
shorter  observation  of  29  days.  A  14  day  time  series  is  adequate  to  provide  the  resol- 
ution of  the  principal  diurnal  constituents,  Kx  and  0„  and  a  15  day  time  series  will  re- 
solve the  principal  semidiurnal  constituents  M2  and  S2.  Some  of  the  longest  tidal 
constituents  require  a  full  19  year  record  for  accurate  resolution  (Schureman,  1958). 

The  current  records  collected  off  Point  Sur  vary  in  length  and  are  not  ideally 
suited  for  the  method  of  linear  combinations.  The  traditional  method  could  provide  only 
a  few  independent  estimates  of  the  tidal  constituents  for  each  month  of  recorded  current 
observations.  Complex  demodulation,  a  least  squares  analysis  technique,  was  used  in- 
stead to  accomplish  the  harmonic  analysis.    The  primary  advantage  realized  with  this 
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Figure  6.  Current  Meter  Record  Filtering:  Time  series  from  mooring  P2  at  100 
m  showing  velocities  in  cm  sec-'  with  LP  filtering,  LLP  filtering  and  the 
band  passed  residual  for  the  u  (top  three  panels)  and  v  (lower  three 
panels)  components. 
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method  is  that  records  of  any  length  can  be  analyzed.  Complex  demodulation  utilizes 
a  set  of  normal  equations  which  are  solved  for  unknown  amplitudes  and  phase  corre- 
sponding to  n  frequency  constituents.  The  number  of  unknowns  solved  for  in  the  set 
of  equations  is  2/7  +  1.  The  number  of  equations  used  in  each  set  is  determined  by  the 
number  of  observations,  m,  utilized  for  each  estimate  and  the  solutions  are  determined 
using  matrix  algebra.   Schematically,  this  process  an  be  represented  as 

[observed  velocities]  =  [known  tidal  constants][unknown  empirical  constants] 

The  least  squares  method  minimizes  the  difference  between  an  observed  current  velocity 
and  the  sum  of  the  solutions  corresponding  to  that  observation.  The  difference,  or  error, 
can  be  attributed  to  unresolved  tidal  constituents,  inertial  currents  or  noise  in  the  signal. 
The  extracted  solutions  reflect  only  those  current  variations  with  a  coherent  phase  and 
no  assumptions  are  made  about  data  outside  the  interval  upon  which  the  solutions  are 
fit  (Pugh,  1986). 

The  analysis  of  diurnal  and  semidiurnal  constituents  was  conducted  starting 
every  six  hours  over  48  hour  segments  and  produced  a  large  number  of  independent  es- 
timates for  each  current  record.  The  48  hour  period  did  not  allow  adequate  separation 
of  the  O,  and  A',  diurnal  and  the  M2  and  S2  semidiurnal  constituents  which  had  prevented 
by  synodic  periods  of  13.661  and  14.765  days  respectively.  Therefore  the  major  tidal 
constituents  (0„  Ku  P{,  Qx,  A/,,  ./,;  M2,  S2,  N2,  and  K2)  were  generalized  as  diurnal  ( 
a)  ~  1  cycle  per  day)  or  semidiurnal  (co  ~  2  cpd).  Longer  period  constituents,  such  as  the 
Lunar  Long  Period  Mm  and  MSf,  were  removed  by  low  pass  filtering.  Harmonics  such 
as  the  Lunar  Terdiurnal  Af3  or  Solar  Terdiurnal  S3  were  not  considered  due  to  the  rela- 
tively small  amplitude  of  these  constituents  when  compared  to  the  primary  diurnal  and 
semidiurnal  constituents.   Higher  order  harmonics  were  removed  by  the  initial  band  pass 
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filtering.  Inclusion  of  the  inertial  frequency  (w  =  1.1  cpd)  as  a  known  frequency  pro- 
vided a  better  fit  of  the  least  squares  solutions  and  reduced  error. 

The  complex  demodulation  program  calculated  the  coefficients  using  a  least 
squares  approximation  to  the  set  of  m  equations  with 

u  or  v  =  Ax  +  A2  sin  o>(/  -  ts)  +  A3  cos  co{i  -  ts)...  (17) 

where  a>  is  the  angular  frequency  of  the  tidal  constituent,  t  is  the  time  of  each  hourly 
record  and  /,  is  the  time  of  the  first  data  point  of  an  entire  current  meter  record.  The 
number  of  terms  can  be  expanded  to  A2n^  to  include  all  desired  constituents  n.  For  u 
velocities,  data  was  placed  in  a  single  row  matrix  U  .  The  coefficients  of  equation  (17) 
are  placed  in  matrix  X  .  The  following  relationship  was  used  to  obtain  the  solutions 
to  the  set  of  equations  represented  by  equation  (17): 


{UxX~l)x{X\X~])  =  AU 


Similar  calculations  were  done  for  the  v  velocities  to  create  matrix  A  V.  Once  the  coef- 
ficients A  were  determined  (the  values  in  matrices  AU  or  AV),  the  amplitude  and  phase 
of  each  frequency  component  was  calculated  using 


Amp„  =  x/(^2J2-(-(/i2fl+1)2 


a  .       f    A2n+\ 

8  =  arctan 


V    2"  / 

With  these  values,  the  current  velocities  can  now  be  defined  as 

uorv=  Amp„  sin(o>(/  —  ts)  +  6) 
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Diurnal  and  semidiurnal  amplitude  and  phase  were  extracted  from  each  current 
meter  record  using  the  complex  demodulation  technique.  An  example  is  provided  in 
Figure  7. 

2.     Construction  of  Current  Ellipses. 

The  amplitude  and  phase  of  clockwise  (QCt  8C)  and  anticlockwise  (Qac,  8ac)  ro- 
tating vectors  were  computed  using  monthly  averages  of  the  component  amplitude  and 
phase  vectors  (U,  8  J  and  (V,8,)  obtained  from  complex  demodulation  using  standard 
equations  (Pugh,  1986): 

Qc  =  \  [U2  +V2-2UV  sin(0v  -  8u)]i 


Qac  =  J  [U2  +  V2  +  2UV  sin(0v  -  8U)]  i 

f  Us'm8u+  Vcos8v  ) 

8C  =  arctan<  — 7. . .  .    a    } 

[  I  cos  8U  —  V  sin  #v  J 

f    (/sinaM+  Vcos8v    } 

8ac  =  arctan<  — .  > 

ac  [  —  I  cos  8U  —  V  sin  8Y  J 


Using  the  rotary  components,  the  dimensions  of  a  current  ellipse  may  be  defined  as 

Semi-major  Axis  =  Qc  +  Qac 

Semi-minor  Axis  =  \QC  —  Qac\ 


Orientation  of  Semi-major  Axis  (clockwise  from  North)  =  -^ r-  (8ac  +  8C) 


Current  rotation  is  clockwise  if  Qc  >  Qac  and  anticlockwise  if  QBe  >  Qc.    Sample  ellipses 
are  shown  for  mooring  P2  during  December  1989  (Figure  8). 
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Figure 


7.      Complex  Demodulation  of  the  Diurnal  Frequency  at  500  m  on  Mooring 

P2:  The  first  two  time  series  plotted  are  amplitudes  in  cm  sec-1  and  the 
next  two  curves  are  the  phases  of  the  u  and  v  diurnal  velocity  compo- 
nents in  degrees  relative  to  the  beginning  of  each  record. 
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TIDAL  CURRENT  ELLIPSE 
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Figure  8.  Tidal  Current  Ellipses:  Current  ellipses  are  shown  for  diurnal  and 
semidiurnal  tidal  currents  at  P2  for  December  1989.  North  is  towards 
the  top  of  the  page. 
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Current  ellipses  from  different  moorings  and  varying  depths  were  evaluated  by 
visual  comparison  of  the  ellipse  orientation,  size,  eccentricity  and  rotation  and  by 
quantitative  comparison  of  the  energy  represented  by  each  ellipse.  Energy  was  com- 
puted in  units  of  cm2  sec-2  by  the  addition  of  the  squared  dimensions  of  the  semi-major 
and  semi-minor  axis. 

C.     SPECTRAL  ENERGY  ANALYSIS 

To  evaluate  the  energy  present  in  the  tidal  currents,  a  program  which  computes  the 
spectral  characteristics  of  a  time  series  was  employed.  The  program,  VCSPC3,  was  used 
to  calculate  the  kinetic  energy  (KE)  spectra  as  the  total  KE  spectra  for  both  the  u  and 
v  velocities  of  a  current  meter  record  over  one  month  periods.  The  program  utilized  a 
Cooley-Tukey  fast  fourier  transform  (FFT).  One  month  records  were  divided  into  five 
sequential,  equal  length  segments.  A  separate  spectrum  was  computed  for  each  segment 
and  the  five  segments  were  then  averaged  to  produce  the  total  energy  spectrum.  A 
parabolic  data  window  was  used  in  the  program  to  eliminate  discontinuities  at  the  be- 
ginning and  end  of  records  created  because  record  lengths  were  not  exact  multiples  of 
the  periods  of  each  wave  present  in  the  signal  {Bendat  and  Pier  sol.  1986).  The  effect  of 
the  data  window  was  to  improve  the  resolution  of  the  narrow  spectral  peaks  in  the  tidal 
frequencies  by  reducing  spectral  side  lobe  leakage  and  associated  distortion.  No  addi- 
tional smoothing  of  the  spectral  estimates  was  done.  The  program  calculated  95% 
confidence  intervals  using  a  Chi  squared  distribution. 

Program  VCSPC3  provided  both  graphical  output  and  data  listings.  The  plots  show 
kinetic  energy  density  in  units  of  cm2  sec-2  cpd-1  versus  frequency  in  cpd.  Confidence 
intervals  are  indicated  by  the  separation  of  parallel  horizontal  lines  above  the  horizontal 
axis.  Distinct  peaks  were  clearly  seen  at  both  the  diurnal  and  semidiurnal  frequencies 
of  these  plots  (Figure  9). 
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Figure  9.  Spectral  Energy  at  P2:  A  representative  plot  from  the  spectral  analysis 
program  is  shown  for  a  one  month  interval.  The  95%  confidence  inter- 
vals is  plotted  below  the  curve. 
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The  listings  show  spectral  density  at  the  diurnal  and  semidiurnal  frequencies  with  a 
frequency  resolution  of  0.05  cpd  and  the  associated  95  percent  confidence  intervals. 

D.     HYDROGRAPHIC  DATA 

A  series  of  temperature  (T),  salinity  (5),  density  (a,),  geostrophic  height  (<t>), 
geostrophic  velocity  (V),  and  spiciness  (n),  values  were  plotted  in  contoured  vertical 
sections  which  passed  through  the  locations  of  moorings  P2  and  P3  (figures  not  shown). 
A  level  of  no  motion  for  the  geostrophic  velocity  calculations  was  chosen  at  1000  m 
except  where  bottom  depths  were  less  than  1000  m  when  the  greatest  common  depth 
between  stations  was  used.  Brunt-Vaisala  (or  buoyancy)  frequency  N  was  computed  for 
each  cruise  at  the  stations  nearest  moorings  P2  and  P3.  Usually  the  CTD  cast  was 
within  a  few  kilometers  of  the  moorings'  assigned  position.  Values  for  Ar  were  deter- 
mined by 


A'2=S 


1    dp 


cz 


(18) 


Values  of  p  used  in  equation  (18)  were  computed  using  the  International  Equation  of 
State  1980  for  seawater  (Miller o  and  Poisson,  1981)  using  S,  T  and  P  values  from  the 
CTD  cast  and  were  smoothed  by  center  averaging  over  10  m  increments.  Calculations 
of  A'2  used  Az  =  10  m.  The  results  from  each  CTD  cast  were  plotted  as  a  vertical  profile 
of  A'2  against  depth.  Most  Ar2  profiles  revealed  the  main  pycnocline  at  depths  between 
50  m  and  200  m.  The  highest  value  of  A',  2.147  x  10_3cps  (a  7.76  minute  period),  was 
found  near  the  main  pycnocline.  Minimum  AT  values  were  found  at  the  bottom  of  each 
cast  with  values  near  2.344  x  lO^cps  (a  71.01  minute  period)  and  were  assumed  to  be 
constant  from  the  deepest  calculated  value  to  the  sea  floor. 

Profiles  of  the  maximum  trapping  frequency  a>c  =  A'sina  were  constructed  with  the 
trapping  frequency  plotted  in  units  of  cpd  against  depth  in  meters.    The  possibility  of 
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bottom  trapping  at  any  particular  CTD  station  is  best  indicated  by  the  magnitude  of  the 
maximum  trapping  frequency  calculated  over  the  lowest  10  m  of  CTD  data  (Figure  10). 
The  deepest  calculated  coc,  typically  50  m  off  the  bottom,  was  extrapolated  to  the  bottom 
using  the  assumption  that  N2  remained  constant.  Any  wave  with  a  frequency  less  than 
or  equal  to  coc  may  be  trapped.  When  coc  is  greater  than  2.0  cpd,  trapping  may  occur 
at  both  diurnal  and  semidiurnal  frequencies.  When  wc  is  less  than  1.0  cpd,  no  trapping 
of  tidal  currents  should  be  expected. 
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Figure   10.         Brunt-  Vaisala     Frequency     and     Trapping     Frequency     Profiles: 

N2  and  A*  sin  a  are  plotted  versus  depth  using  hydrographic  data  ob- 
tained near  P3  in  March  1990. 
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IV.     RESULTS 

A.     CURRENT  ELLIPSES 

Monthly  averaged  tidal  current  ellipses  were  constructed  using  the  amplitudes  and 
phases  obtained  from  the  complex  demodulation  of  current  meter  records  (Figures  11 
to  14).  The  plotted  ellipses  display  a  great  deal  of  variability  in  their  size,  orientation 
and  direction  of  rotation.  The  size  and  energy  represented  by  the  current  ellipses  was 
found  to  be  greatest  for  both  tidal  frequencies  at  the  deepest  current  meter  of  mooring 
P2. 

1.     Diurnal  Current  Ellipses 

Diurnal  current  ellipses  for  mooring  P2  were  constructed  for  seventeen  consec- 
utive months  (Figures  11  and  12).  The  energy  represented  by  the  ellipses  was  calculated 
using  the  sum  of  the  squared  axis  lengths  to  compute  kinetic  energy  in  units  of 
cm2  sec-2  .  At  the  100  m  level  the  average  energy  of  the  ellipses  was  18.0  cm2  sec-2.  The 
direction  of  rotation  of  76%  of  these  ellipses  was  counterclockwise.  At  350  m,  the  av- 
erage energy  dropped  to  14.5  cm2  sec-2  with  70%  of  the  ellipses  rotating  counterclock- 
wise. At  500  m,  the  average  ellipse  energy  increased  to  22.3  cm2  sec-2  with  65%  rotating 
counterclockwise.  The  diurnal  current  ellipses  were  larger  and  more  energetic  at  500  m 
than  at  350  m  for  each  of  the  seventeen  months  analyzed  and  there  were  twelve  months 
when  more  energy  was  present  at  500  m  than  at  100  m.  During  July  1989  the  diurnal 
current  rotation  at  all  three  P2  current  meters  was  clockwise  and  there  were  eight 
months  when  the  rotation  at  three  depths  was  uniformly  counterclockwise. 

The  current  meter  records  for  six  months  were  used  to  construct  diurnal  current 
ellipses  at  mooring  P3  (Figure  12).  At  the  100  m  level  the  average  energy  represented 
by  the  ellipses  was  25.1  cm2  sec-2  with  66%  of  the  ellipses  rotating  counterclockwise. 
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Figure   11.      Diurnal  Current  Ellipses  from  May  1989  to  November  1989  at  Mooring 

P2:  Horizontal  ellipses  are  shown  with  north  up  and  the  direction  of 
rotation  indicated  by  arrows.  A  velocity  scale  is  indicated  at  the  bot- 
tom. 
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Figure  12.      Diurnal  Current  Ellipses  from  December  1989  to  September  1990  at 
Moorings  P2  and  P3:     Details  are  as  in  Figure  11. 
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Figure  13.       Semidiurnal  Current  Ellipses  from  May   1989  to  November  1989  at 
Moorings  P2  and  P3:     Details  are  as  in  Figure  11. 
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Semidiurnal  Current  Ellipses  from  December  1989  to  September  1990 
at  Mooring  P2  and  P3:     Details  are  as  in  Figure  1 1. 
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At  350  m,  the  average  energy  dropped  to  12.4  cm2  sec-2  .  At  500  m,  four  months  were 
analyzed  and  resulted  in  an  average  energy  of  13.5  cm2  sec-2  with  75%  rotating 
counterclockwise.  The  smallest  ellipses  and  lowest  energy  (9.5  cm2  sec-2)  were  found  at 
1000  m  with  66%  of  the  ellipses  rotating  counterclockwise.  The  ellipses  were  larger  and 
contained  more  energy  at  500  m  than  did  those  at  either  350  m  or  1000  m.  There  were 
four  months  when  the  current  rotation  at  each  P3  current  meter  was  counterclockwise. 

The  bathymetric  contours  at  mooring  P2  are  oriented  at  approximately  005  °T 
and  the  average  orientation  of  the  500  m  diurnal  tidal  current  ellipses  at  P2  was  049° 
T.  The  semi-major  axis  of  the  diurnal  ellipses  crossed  the  isobaths  a  mean  relative  angle 
of  44°.  At  mooring  P3  the  bathymetry  is  oriented  at  020°  T.  The  average  orientation 
of  diurnal  tidal  ellipses  at  1000  m  was  074°  T  and  the  semi-major  axis  of  these  ellipses 
crossed  the  isobaths  at  a  mean  relative  angle  of  54°. 
2.     Semidiurnal  Current  Ellipses. 

The  semidiurnal  current  ellipses  at  P2  and  P3  were  larger  and  more  energetic 
than  those  at  diurnal  frequencies.  The  orientation  of  the  semidiurnal  ellipses  was  less 
variable  than  the  diurnal  ellipses  with  almost  all  of  them  oriented  in  a  north-northeast 
direction. 

At  mooring  P2,  the  seventeen  month  average  energy  represented  by  the  ellipses 
at  100  m  was  59.1  cm2  sec-2.  The  rotation  of  58%  of  these  ellipses  was  counterclockwise. 
At  350  m,  the  average  energy  decreased  to  48.2  cm2  sec-2  and  52%  of  the  currents  ro- 
tated counterclockwise.  At  500  m  the  average  energy  increased  to  85.4  cm2  sec-2  with 
47%  of  the  ellipses  displaying  counterclockwise  rotation.  The  semidiurnal  currents  dis- 
played the  same  trends  exhibited  by  the  diurnal  current  ellipses  where  the  largest  ellipses 
and  greatest  energies  were  found  at  the  deepest  level  and  the  smallest  ellipses  and  least 
energy  was  at  350  m.    The  rotation  of  semidiurnal  current  ellipses  at  P2  was  clockwise 
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at  all  three  current  meters  during  four  months  and  counterclockwise  during  four  other 
months.   The  remaining  nine  months  contained  ellipses  with  opposing  rotation. 

The  semidiurnal  current  ellipses  at  mooring  P3  were  notably  smaller  than  those 
at  P2.  The  average  energy  represented  by  the  ellipses  at  100  m  was  28.4  cm2  sec-2  with 
66%  rotating  clockwise.  At  350  m  the  average  energy  decreased  to  22.7  cm2  sec-2  with 
66%  of  the  currents  rotating  counterclockwise.  A  four  month  average  at  500  m  resulted 
in  an  average  semidiurnal  energy  of  34.5  cm2  sec-2  with  all  four  ellipses  rotating  clock- 
wise. At  1000  m  the  average  energy  was  only  8.3  cm2  sec-2.  At  this  level,  66%  of  the 
ellipses  rotated  clockwise.  Like  the  diurnal  currents  at  P3,  the  semidiurnal  currents  were 
greatest  at  500  m.  The  smallest  ellipses  and  least  energy  was  at  the  1000  m  level.  At 
mooring  P3  two  of  the  six  months  analyzed  had  semidiurnal  currents  at  four  depths  ro- 
tating in  a  clockwise  direction  and  there  was  one  month  when  the  rotation  was 
counterclockwise  at  each  depth. 

The  average  orientation  of  500  m  semidiurnal  current  ellipses  at  P2  was  030°  T 
and  the  isobaths  were  crossed  at  an  angle  of  25°.  The  average  orientation  of  1000  m 
semidiurnal  current  ellipses  at  P3  was  043°  T  and  the  semi-major  axis  crossed  the 
isobaths  at  23°. 

B.     ENERGY  PROFILES 

Vertical  profiles  of  energy  at  both  the  diurnal  and  semidiurnal  frequencies  (Figures 
15  and  16)  were  constructed  for  each  month  using  the  output  from  the  spectral  analysis 
program.  These  plots  contain  information  similar  to  that  of  the  tidal  current  ellipses 
expressed  in  a  different  format.  At  P2  the  highest  energy  levels  were  at  the  deepest 
current  meter  for  both  frequencies.  The  energy  levels  at  P3  show  low  energy  levels  at 
the  deeper  current  meters. 
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Figure  15.      Spectral  Energy  Profiles  from  May  1989  to  November  1989  at  Mooring 

P2:     All  plots  are  scaled  identically. 
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Figure  16.      Spectral  Energy  Profiles  from  December  1989  to  September  1990  at 
Moorings  P2  and  P3:     All  plots  are  scaled  identically. 
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1.  Diurnal  Energy  Profiles 

At  mooring  P2  the  average  monthly  diurnal  energy  levels  were  13.9  cm2  sec-2  at 
100  m,  13.4  cm2  sec-2  at  350  m  and  20.2  cm2  sec2  at  500  m.  The  energy  at  500  m  ex- 
ceeded that  at  350  m  for  each  of  the  seventeen  months  analyzed  except  September  1989. 
The  energy  at  500  m  was  also  higher  than  that  at  100  m  except  for  the  months  of  June 
1989,  February  1990  and  April  1990.  The  curves  plotted  in  Figures  15  and  16  exhibit  a 
C  shape  with  minimum  energy  levels  typically  found  at  a  depth  of  350  m. 

Diurnal  energy  at  P3  was  greatest  near  the  surface.  Average  monthly  values 
were  17.3  cm2  sec-2  at  100  m.  7.2  cm2  sec-2  at  350  m  9.7  cm2  sec-2  at  500  m  and  9.8 
cm2  sec-2  at  1000  m.  The  curves  plotted  below  100  m  are  nearly  uniform  and  show  vir- 
tually no  amplification  of  energy  toward  the  bottom. 

2.  Semidiurnal  Energy  Profiles 

The  semidiurnal  current  energy  levels  at  P2  were  likewise  greatest  at  the  deepest 
current  meters.  The  same  'C  shaped  curve  observed  with  diurnal  currents,  with  a  min- 
imum energy  at  350  m,  was  exhibited  by  the  semidiurnal  currents.  Average  monthly 
values  of  63.1,  47.6,  and  77.4  cm2  sec2  were  found  at  100  m,  350  m  and  500  m  respec- 
tively. The  energy  at  500  m  was  greater  than  that  at  350  m  for  fourteen  of  the  seventeen 
months  analyzed.  The  energy  at  500  m  exceeded  that  at  100  m  for  nine  of  the  seventeen 
months. 

At  P3,  the  average  semidiurnal  energy  levels  were  fairly  consistent  until  drop- 
ping off  sharply  at  1000  m.  Values  of  20.1,  21.3,  30.3  and  7.3  cm2  sec-2  were  found  at 
100  m,  350  m,  500  m  and  1000  m.  No  bottom  amplification  of  semidiurnal  energy  was 
evident  at  P3. 

C.     TEMPORAL  VARIABILITY 

The  same  data  used  to  plot  the  vertical  profiles  were  used  to  create  time  series 
showing  the  change  of  energy  at  each  depth  for  the  two  moorings.    Diurnal  bottom 
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trapping  was  evident  at  P2,  with  energy  at  500  m  greater  than  the  energy  at  350  m  for 
the  entire  seventeen  month  period  analyzed  (Figure  17).  Bottom  trapping  was  always 
present  and  was  most  pronounced  during  July  -  August  1989,  October  1989  -  January 
1990  and  May  -  September  1990  when  the  energy  at  500  m  exceeded  that  of  any  other 
depth.  The  magnitude  of  the  diurnal  energy  at  500  m  was  statistically  greater  than  350 
m  or  100  m  at  95  percent  confidence  for  the  period  October  1989  through  January  1990. 
At  mooring  P3  no  diurnal  bottom  trapping  was  indicated  by  the  energy  at  350  m,  500 
m  or  1000  m  (Figure  18).  Diurnal  energy  did  not  increase  toward  the  bottom  and  the 
energy  near  the  surface,  at  100  m,  was  always  greater  than  at  deeper  current  meters. 

Bottom  trapping  of  semidiurnal  energy  occurred  at  P2  for  fourteen  of  seventeen 
months  when  the  energy  at  500  m  exceeded  the  energy  at  350  m  (Figure  19).  In  addition, 
the  energy  at  500  m  was  greater  than  100m  for  nine  of  the  seventeen  months  analyzed. 
Bottom  trapping  at  500  m  was  statistically  significant  using  a  95  percent  confidence  in- 
terval for  September  1989,  December  1989,  April  -  May  1990  and  September  1990.  At 
P3  no  semidiurnal  trapping  was  evident  and  the  energy  at  1000  m  was  less  than  at  any 
other  depth  for  five  of  the  six  months  analyzed  (Figure  20). 

D.     APPLICATION  OF  TRAPPING  THEORY 

The  size  of  the  current  ellipses  and  energy  distribution  at  each  depth  suggest  that 
an  amplification  of  energy  occurs  near  the  bottom  of  mooring  P2.  A  similar  pattern  was 
not  obvious  at  P3.  To  compare  observed  motions  to  theoretical  trapped  waves  an  ex- 
ponential curve,  representing  the  decay  of  trapped  energy,  was  compared  to  the  observed 
profile  of  kinetic  energy.  The  theoretical  curves  were  constructed  using  equation  (12) 
with  k  given  by  equation  (16). 

1.     Trapping  Scales  and  Decay  Curves 

At  mooring  P2  the  maximum  trapping  frequency  is  exactly  equal  to  the  diurnal 
frequency  (coc  =  1  cpd)  when,  by  equation  (7),  N  =  1.930  x  10-4  cps  .    Using  H  =  800  m 
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Time  Series  of  Diurnal  Energy  at  Mooring  P2:  The  temporal  vari- 
ability of  kinetic  energy  is  shown  with  curves  representing  the  energy 
levels  at  current  meter  depths. 
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Figure  18.      Time  Series  of  Diurnal  Energy  at  Mooring  P3: 
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Figure  17. 
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and   /=  1.375  x  10-5  cps    in    equation    (16)    results    in    a    trapping    coefficient    of 

k  =  3.3116  x  10~3  nr1.   The  horizontal  length  scale  of  this  wave  is  L  =  —r-=  26.632  km. 

k 

Similar  calculations  were  conducted  for  the  first  five  baroclinic  modes  at  both  P2  and 
P3  for  diurnal  and  semidiurnal  trapping  using  N  values  determined  by  equation  (7).  The 
coefficients  are  tabulated  in  Tables  3  and  4. 


Table  3.     TRAPPING  COEFFICIENTS  AT  MOORING  P2 


v)c  =  1  cpd 

w(  =  2  cpd 

k,  =  3.311xlO-3m-' 

k,  =  1.472xl0-3m-' 

K2  =  6.622xl0-3nr1 

K:  =  2.943xl0-3m-' 

K3  =  9.934xl0-3m-1 

K3  =  4.415xl0-3m-' 

k4=  1.321x10-*  nr1 

K4  =  5.887xl0-3m-' 

ks=  1 .  652x1 0-2m-' 

K'5  =  7.358xlO-3m-' 

Table  4.     TRAPPING  COEFFICIENTS  AT  MOORING  P3 


vjc  =  1  cpd 

v)c  =  2  cpd 

k}  =6.622xl0-3  m-1 

K,  =  2.949xl0-3m-1 

k2=  1.324xl0-2m-' 

K2  =  5.887xlO-3m-' 

k3  =  1.986xl0-2m-] 

k3  =  8.830x10-3  m-1 

k4  =  2.643xl0-2  m-1 

k4=  1.177xl0-3m-1 

k<  =  3.33 lxlO-2  nr1 

K{=  1.476xl0-2  m-1 

For  the  baroclinic  modes,  the  trapping  coefficients  in  Tables  3  and  4  were  used 
to  produce  a  vertical  trapping  scale.  The  trapping  scale  (or  e-folding  scale)  indicates  the 
height  off  the  bottom  where  the  energy  of  a  trapped  wave  decayed  by  a  factor  of  e~x  . 
Most  of  the  trapped  energy  for  each  wave  was  contained  in  a  bottom  layer  with  a 
thickness  (determined  by  the  trapping  scale)  indicated  in  Table  5.  Note  that  the  two 
deepest  current  meters  at  P2  (located  450  m  and  300  m  from  the  bottom)  fell  within  this 
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layer  but  the  lowest  instrument  at  P3  (located  800  m  from  the  bottom)  did  not.  For 
higher  modes  trapping  is  stronger,  trapping  scales  are  smaller  and  energy  decays  more 
rapidly  from  the  bottom. 

Table  5.     TRAPPING  SCALES  AT  MOORINGS  P2  AND  P3 


Mooring 

wc  =  1  cpd 

v>c  —  2  cpd 

P2 

302  m 

679  m 

P3 

151  m 

339  m 

Using  the  appropriate  trapping  coefficient,  exponential  decay  curves  for  each 
mooring  were  constructed  which  represent  the  first  two  baroclinic  modes  at  diurnal  and 
semidiurnal  frequencies  at  moorings  P2  and  P3  (Figures  21  and  22).  Similar  curves  were 
fit  to  observed  spectral  energy  profiles  to  determine  if  the  vertical  energy  distribution  off 
Point  Sur  indicated  the  presence  of  bottom  trapped  waves.  The  curves  were  constructed 
using  A'  values  determined  from  CTD  data  whenever  the  diurnal  or  semidiurnal  fre- 
quencies were  less  than  the  maximum  trapping  frequency. 
2.     Analysis  of  Bottom  Trapping 

CTD  data  from  each  of  the  hydrographic  cruises  in  Table  2  was  analyzed  to  see 
if  bottom  trapping  was  possible  at  P2  or  P3  at  diurnal  or  semidiurnal  frequencies  using 
the  criteria  that  ojc  =  AT  sin  a  must  be  equal  to  or  greater  than  the  tidal  frequencies. 
Each  CTD  cast  taken  at  a  mooring  site  was  assumed  to  be  representative  of  the  aver- 
aging interval  used  to  compute  the  spectral  estimates  (one  month). 

When  bottom  trapping  was  possible,  theoretical  decay  curves  were  constructed 
using  observed  values  of  H,  N,  and  appropriate  trapping  coefficients  (k).  The  curves 
were  compared  to  the  observed  energy  levels  obtained  from  spectral  analysis  by  match- 
ing the  curves  at  the  level  of  the  deepest  current  meter  (500  m  at  P2)  and  comparing  the 
theoretical  and  observed  energy  at  other  (100  m  and  350  m  at  P2)  levels.    When  the 
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Figure  21.  Theoretical  Trapping  Profiles  at  Mooring  P2:  The  exponential  decay 
of  trapped  energy  (with  arbitrary  units)  is  illustrated  for  the  first  and 
second  modes  at  tidal  frequencies. 

modeled  energy  fell  within  the  95  percent  confidence  intervals  of  the  observations,  the 

change  of  energy  with  depth  was  considered  to  be  consistent  with  bottom  trapped  waves. 

The  buoyancy  frequencies  AT,  maximum  trapping  frequencies  coe  and  trapping 

coefficients  k,  for  each  month  when  data  were  available  are  tabulated  in  Table  6.    The 
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Figure  22.  Theoretical  Trapping  Profiles  at  Mooring  P3:  The  exponential  decay 
of  trapped  energy  (with  arbitrary  units)  is  illustrated  for  the  first  and 
second  modes  at  tidal  frequencies. 

maximum  trapping  frequency  coc  at  P2  exceeded  the  diurnal  frequency  (co  =  1.0  cpd) 
each  of  the  eight  months.  The  maximum  trapping  frequency  exceeded  the  semidiurnal 
frequency  (w  =  2.0  cpd)  except  for  the  month  of  January  1990  (w£  =  1.22  cpd). 
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Table  6.     OBSERVED  BUOYANCY  FREQUENCIES,  TRAPPING  FREQUEN- 
CIES AND  MODE  1  TRAPPING  COEFFICIENTS  AT  MOORING  P2 


Month 

A' (cpd) 

o)c  (cpd) 

KDitmul    (m       ) 

KStmidhrnal   \V^      ) 

MAY  89 

66.2 

3.97 

3.30555  x  10~3 

6.611 16  x  10-3 

JUL  89 

43.0 

2.57 

3.30556  x  10-3 

6.61127  x  10-3 

SEP  89 

39.4 

2.36 

3.30557  x  10-3 

6.61130  x  10-3 

NOV  89 

34.7 

2.08 

3.30558  x  10-3 

6.61137  x  10-3 

JAN  90 

20.3 

1.22 

3.30564  x  10-3 

No  Trapping 

MAR  90 

51.8 

3.10 

3.30556  x  10-3 

6.61121  x  10-3 

MAY  90 

38.7 

2.32 

3.30557  x  10-3 

6.61131  x  10-3 

J  UN  90 

55.7 

3.34 

3.30556  x  10-3 

6.61120  x  10-3 

AVERAGE 

43.7 

2.62 

3.30557  x  10-3 

6.61 126  x  10-3 

The  diurnal  energy  distribution  at  P2  (Figures  23  and  24)  agreed  with  the  mod- 
eled decay  for  seven  of  the  eight  months.  During  January  1990,  the  observed  energy  at 
350  m  was  less  than  indicated  by  theory.  The  theoretical  curves,  extrapolated  to  a  depth 
of  800  m,  indicated  an  average  kinetic  energy  level  of  40.5  cm2  sec-2  trapped  at  the  sea 
floor. 

The  maximum  trapping  frequency  oc  at  P2  exceeded  the  semidiurnal  frequency 
(cue=  2.0  cpd)  during  seven  of  the  eight  months.  The  observed  energy  at  100  m  during 
May  1989  was  higher  than  was  modeled  by  matching  the  theoretical  curve  to  the  500 
m  observation.  The  modeled  curve  fell  within  the  95  percent  confidence  intervals  for  the 
six  remaining  months.  Extrapolation  of  the  theoretical  curve  to  800  meters  produced 
an  average  kinetic  energy  level  of  310.0  cm2  sec-2  at  the  bottom  over  the  time  period 
studied. 

The  buoyancy  frequency  calculated  over  the  deepest  10  m  of  CTD  data  ranged 
from  20.3  cpd  to  66.2  cpd  and  the  magnitude  of  the  maximum  trapping  frequency,  which 
was  proportional  to  N  ,  ranged  from  1.22  cpd  to  3.97  cpd.  The  fluctuations  of  A7  had  a 
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much  smaller  effect  on  the  trapping  coefficients.  The  diurnal  trapping  coefficient  at  P2 
varied  only  9.0  x  10~8  m_1.  The  semidiurnal  trapping  coefficient  had  a  range  of  only 
2.1  x  lO-7  nr1.  The  consistent  trapping  coefficients  produced  nearly  identical  trapping 
scales  at  P2.  The  diurnal  trapping  scale  (defined  by  the  e-folding  scale)  was  302  m  and 
the  semidiurnal  trapping  scale  was  151  m  for  each  of  the  eight  months  analyzed. 

The  diurnal  energy  distribution  at  P3  did  not  have  a  profile  indicative  of  bottom 
trapped  waves.  Since  the  maximum  trapping  scale  for  a  diurnal  wave  at  P3  was  ap- 
proximately 151  m  and  the  deepest  current  meter  at  P3  was  800  m  above  the  sea  floor, 
any  energy  from  a  trapped  wave  would  have  decayed  99.5  percent  before  reaching  this 
current  meter  and  could  not  be  observed.  The  maximum  trapping  frequency  at  this 
mooring  never  exceeded  cdc  —  2.0  cpd  and  no  semidiurnal  trapping  occurred. 
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Figure  23. 


Trapping  Frequencies  and  Energy  Profiles  at  P2  During  1989:  Maxi- 
mum trapping  frequencies  are  plotted  against  depth  with  the  dashed 
line  indicating  2.0  cpd  in  the  top  row.  The  middle  row  shows  observed 
diurnal  spectral  energy  profiles  (solid  lines)  with  95%  confidence  inter- 
vals (horizontal  bars).  Theoretical  energy  profiles  (dashed  lines),  based 
on  a  decay  scale  calculated  using  .V  from  the  bottom  ofCTD  casts  (or 
at  800  m  for  July  1989)  ,  are  matched  to  the  observations  at  500  m. 
Similar  profiles  are  shown  in  the  third  row  for  semidiurnal  energy. 
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Figure  24.      Trapping  Frequencies  and  Energy  Profiles  at  P2  During  1990:     Data  is 

plotted  as  described  in  Figure  23.  No  theoretical  semidiurnal  energy 
profile  is  shown  for  January  1990  when  A*  sin  a  near  the  bottom  was 
less  than  2.0  cpd. 
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V.     DISCUSSION 

A.     BOTTOM  TRAPPED  WAVES 

The  trapping  or  amplification  of  tidal  band  kinetic  energy  near  the  ocean  bottom 
was  evident  from  the  data  obtained  at  mooring  P2.  Both  the  current  ellipses  and  the 
energy  distribution  displayed  a  significant  increase  from  a  depth  of  350  m  to  500  m 
during  the  seventeen  months  of  data  analyzed.  A  similar  increase  was  not  evident  from 
the  data  obtained  at  mooring  P3. 

The  two  methods  of  analysis  used  to  examine  the  currents  produced  results  which 
were  in  general  agreement.  The  energy  levels  derived  through  complex  demodulation 
and  current  ellipses  were  larger  in  magnitude  than  those  obtained  through  spectral 
analysis.  The  total  KE  represented  by  the  ellipses  was  10.7  %  greater  than  the  energy 
obtained  from  spectral  analysis.  The  lower  energy  spectral  analysis  values  may  be  due 
to  Dopplcr  shifting  of  small  wavelength  internal  tides  to  other  parts  of  the  spectrum. 
This  shifting  could  spread  energy  around  the  tidal  band  so  that  energy  measured  pre- 
cisely at  the  frequency  of  a  tide,  using  the  high  resolution  FFT  of  VCSPC3,  improperly 
excluded  tidal  energy  located  in  sidebands  {IV wish,  1975).  The  difference  between  the 
energy  from  ellipses  and  spectral  analysis  was  greatest  at  P3  during  December  1989.  The 
spectral  density  plots  for  the  current  meter  records  at  100  m  and  350  m  show  that  the 
diurnal  peak  was  spread  out  over  a  wider  frequency  band  (perhaps  by  Doppler  shifting) 
and  had  a  smaller  amplitude  than  for  other  months  of  the  same  record. 

Both  ellipses  and  spectral  energy  profiles  indicated  more  energy  was  present  at  100 
m  than  at  350  m  at  both  moorings  for  diurnal  and  semidiurnal  currents.  The  current 
meters  at  100  m  were  above  the  main  pycnocline  and  the  higher  energy  levels  in  the 
surface  layer  could  be  the  result  of  wind  forcing  at  or  near  tidal  frequencies.    Strong 
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coastal  winds  with  diurnal  variations  have  often  been  observed  along  the  California 
coast  (Beardsley  et  ai,  1987).  The  pycnocline  could  prevent  coupling  of  surface  energy 
to  deeper  layers.  Diurnal  tidal  currents  may  also  be  influenced  by  inertial  currents  with 
a  frequency  is  1.1  cpd  off  Point  Sur.  The  narrow  frequency  separation  between  diurnal 
and  inertial  currents  could  result  in  smearing  or  leakage  of  energy  into  the  diurnal  band. 
The  observed  distribution  of  tidal  energy  had  a  vertical  decay  rate  indicative  of 
bottom  trapping  at  mooring  P2.  The  fit  of  a  theoretical  curve  to  the  data  fell  within  the 
95  percent  confidence  intervals  with  few  exceptions.  In  January  1990  the  observed 
diurnal  energy  at  350  m  was  less  than  predicted  by  bottom  trapped  theory  using  the 
buoyancy  frequency  calculated  from  a  single  CTD  cast.  If  the  buoyancy  frequency  had 
been  greater  the  trapping  coefficient,  k,  would  have  been  larger.  A  larger  k  would  result 
in  a  faster  decay  rate  and  the  theoretical  energy  at  350  m  would  have  been  closer  to 
observations.  The  semidiurnal  energy  at  100  m  during  May  1989  was  higher  than  ex- 
pected using  the  modeled  decay  curve  and  these  energetic  semidiurnal  currents  may  re- 
sult from  surface  (wind)  forcing.  Another  factor  in  the  vertical  semidiurnal  energy 
distribution  is  free  internal  waves  which  exist  at  superinertial  frequencies.  Free  internal 
waves  could  be  generated  by  tidal  forcing  across  the  continental  shelf  break  and  could 
then  propagate  into  deeper  stratified  water  (Baines,  1982).  Horizontal  velocities  over  the 
shelf  have  been  observed  where  the  highest  velocities  were  near  the  surface  and  bottom. 
The  velocity  structure  was  attributed  to  the  first  dynamical  mode  {Rosenfeld,  1990)  and 
the  same  structure  may  dominate  internal  tides  over  the  continental  slope.  The 
semidiurnal  currents  at  P2,  located  close  to  the  shelf  break,  were  twice  as  energetic  as  the 
semidiurnal  currents  found  further  offshore  at  P3  and  did  have  more  energy  at  the  top 
and  bottom  current  meters  (100  m  and  500  m)  than  at  the  middle  (350  m).  Bottom 
trapping  did  not  occur  with  any  discernable  seasonal  pattern.    The  absence  of  distinct 
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temporal  variation  agrees  with  other  observations  of  internal  tides  along  the  California 
coast  (Soble  et  a!.,  1987). 

The  theoretical  curves  plotted  for  P2  in  Figures  23  and  24  were  extrapolated  to  800 
m.  The  average  diurnal  energy  represented  by  these  estimates  was  40.5  cm2  sec-2.  Using 
the  average  diurnal  current  ellipse  eccentricity  (e  =  3.42),  diurnal  currents  at  800  m 
would  have  had  a  semi-major  axis  of  6.1  cm  sec-1  and  a  semi-minor  axis  of 
1.8  cm  sec-1.  The  theoretical  energy  for  semidiurnal  currents  extrapolated  to  800  m  was 
310.0  cm2  see-2.  An  average  semidiurnal  eccentricity  of  e  =  3.26  produced  an  ellipse  re- 
presenting bottom  currents  with  semi-major  and  semi-minor  axis  lengths  of  16.8 
cm  sec-1  and  5.2  cm  sec-1  respectively.  The  magnitude  of  these  currents  is  reasonable, 
however,  bottom  trapped  theory  assumes  inviscid  boundary  conditions  and  bottom 
friction  may  reduce  the  velocity  of  predicted  bottom  trapped  currents. 

The  average  wavelength  of  bottom  trapped  diurnal  waves,  computed  using  the  av- 
erage trapping  coefficient  in  Table  5,  was  66.5  km.  The  average  wavelength  of  bottom 
trapped  semidiurnal  waves  was  33.2  km.  This  compares  favorably  to  first  baroclinic 
mode  internal  waves  of  semidiurnal  frequency  observed  during  CODE  which  had  wave- 
lengths between  20  and  30  km  (Rosenfeld,  1990). 

The  rotation  of  tidal  current  ellipses  is  not  described  by  bottom  trapped  theory.  The 
observed  rotation  of  the  current  ellipses  was  variable  with  73  percent  of  the  diurnal  el- 
lipses rotating  counterclockwise  and  52  percent  of  the  semidiurnal  ellipses  rotating 
clockwise.  Barotropic  Kelvin  waves,  free  internal  waves  and  shelf  waves  all  influence  the 
rotation  of  current  ellipses  at  tidal  frequencies.  The  observed  direction  of  rotation  is 
determined  by  the  cumulative  effect  of  these  unresolved  wave  regimes. 

Barotropic  semidiurnal  tides  are  dominated  by  the  M2  constituent.  The  horizontal 
velocities  generated  by  Kelvin  waves  at  this  frequency  are  normally  aligned  with  topog- 
raphy (Noble  et  ai,  1987).    The  observed  tidal  current  ellipses  all  had  cross-slope  com- 
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ponents  and  vertical  velocity  structures  which  may  reflect  the  influence  of  a  baroclinic 
wave  structure.  At  the  deepest  current  meters,  the  diurnal  current  ellipses  were  oriented 
more  across  the  bottom  slope  (P2  =  44°,  P7>  =  54°)  than  were  the  semidiurnal  ellipses  ( 
P2  =  25°,  P3  =  23°  ).  With  this  geometry  the  particle  motion  of  diurnal  currents  had  a 
larger  upslope  component  and  more  bottom  trapping  than  semidiurnal  currents  of  the 
same  magnitude.  Even  though  observed  diurnal  currents  were  significantly  less  energetic 
than  semidiurnal  currents,  diurnal  bottom  trapping  was  always  evident  at  mooring  P2. 

Both  moorings  had  a  similar  bottom  slope  a  and  the  currents  analyzed  were  limited 
to  diurnal  and  semidiurnal  frequencies.  The  only  variable  in  determining  the  limiting 
trapping  frequency  (  <yc  =  N  sin  a  )  was  the  buoyancy  frequency.  Higher  values  of  N 
resulted  in  higher  trapping  frequencies.  For  each  of  the  eight  months  when  CTD  casts 
were  made  at  P2,  diurnal  trapping  was  possible  as  indicated  by  the  calculated  values  of 
(oe  and  semidiurnal  trapping  was  indicated  during  seven  of  these  months.  The  current 
ellipses  and  energy  profiles  exhibited  bottom  trapping  in  each  case.  No  bottom  trapping 
was  evident  at  P3  although  there  were  two  months  with  o>c  values  high  enough  to  permit 
the  trapping  of  diurnal  energy.  The  absence  of  any  notable  bottom  trapping  at  P3  may 
be  due  to  the  trapping  scale  and  the  depth  of  the  instruments  at  that  mooring.  All  three 
current  meters  at  P3  were  above  the  e-folding  depths  of  bottom  trapped  diurnal  and 
semidiurnal  energy.  If  any  trapped  energy  was  present,  the  magnitude  of  the  signal  was 
indiscernible  at  the  P3  current  meters. 

The  difference  in  maximum  trapping  frequencies  at  each  mooring  is  caused  by  dif- 
ferent hydrographic  conditions  at  P2  and  P3.  The  cruise  of  November  1989  can  be  used 
to  illustrate  typical  conditions.  Density  gradients  are  greater  near  the  ocean  bottom  at 
P2  than  at  P3  (Figure  25)  and  the  gradients  produce  higher  N  values  at  P2  than  those 
calculated  at  P3.  The  stronger  gradients  at  P2  permit  trapping  of  waves  with  frequencies 
higher  than  the  maximum  frequency  of  trapped  waves  at  mooring  P3. 
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Figure  25.       Vertical  Density  Section  Off  Point  Sur,  California  During  November, 
1989:     The  positions  of  moorings  P2  and  P3  are  indicated. 

B.     LIMITATIONS 

Several  assumptions  were  made  in  comparing  the  observations  at  P2  with  theoretical 
predictions  for  bottom  trapped  waves.  The  first  was  that  the  tidal  currents  were  always 
oriented  across  the  sea  floor  slope  so  that  the  angle  </>  =  90°.  This  angle  permits  maxi- 
mum trapping  and  simplifies  the  equation  o  =  N  sin  a  sin  4>  to  coc  =  N  sin  a.  At  angles 
less  than  90°,  the  limiting  trapping  frequency  oc  will  be  lower  and  the  bottom  trapping 
of  tidal  energy  may  not  occur.   If  </>  was  assumed  to  be  the  angle  at  which  the  semi-major 
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axis  of  the  tidal  current  ellipses  crossed  the  bottom  contours,  then  the  maximum  trap- 
ping frequency  a  at  P2  (</>  =  44°  at  500  m)  would  still  exceed  1.0  cpd  for  seven  of  the 
eight  months  when  N  values  were  available.  For  the  P2  semidiurnal  ellipses  ( 
4>  =  25°  at  500  m),  a  would  not  have  exceeded  2.0  cpd  for  any  month.  If  no  semidiurnal 
trapping  was  possible,  then  the  vertical  distribution  of  energy  at  semidiurnal  frequencies 
must  be  attributed  to  other  factors  including  free  internal  waves. 

A  second  assumption  was  that  a  single  CTD  cast  provided  information  which  could 
be  used  to  interpret  monthly  averaged  current  meter  data.  The  CTD  casts  typically 
ended  50  m  above  the  sea  floor  and  salinity  and  temperature  values  over  the  deepest  10 
m  interval  were  used  to  calculate  N  for  each  month.  These  N  values  were  assumed  to 
be  constant  from  the  bottom  of  the  CTD  cast  to  the  sea  floor.  Trapping  frequencies 
calculated  with  N  were  similarly  extrapolated.  The  use  of  extrapolated  values  was  sup- 
ported by  the  CTD  cast  of  July  1989  which  was  slightly  offshore  of  mooring  P2  and  ex- 
tended to  903  m  (the  10  m  interval  centered  at  800  m  was  used  to  calculate  N  at  P2). 
The  values  of  N  were  fairly  constant  during  July  1989  (42.2  cpd  -  40.9cpd)  at  the  depths 
where  extrapolation  was  used  (650  m  to  800  m)  for  the  seven  other  CTD  casts.  As  dis- 
cussed in  the  previous  section,  variations  in  TV  could  affect  the  fit  of  the  theoretical  decay 
curves  to  observations.  The  a>c  and  k  values  using  a  constant  N  were  compared  to 
(oc  and  k  computed  using  Ar  obtained  by  extending  the  slope  of  the  N2  curves  to  800  m. 
The  differences  created  slight  variations  of  coc  which  did  not  conflict  with  any  of  the  re- 
sults presented  in  Chapter  IV.  The  differences  in  the  k  values  were  very  small  and  con- 
sidered insignificant. 

The  instrumentation  at  each  of  the  moorings  was  not  ideally  positioned  for  the  study 
of  bottom  trapped  waves.  The  current  meters  were  too  sparse  to  provide  a  good  profile 
of  the  currents  throughout  the  water  column.  The  current  meters  at  the  100  m  level  were 
above  the  pycnocline  and  may  have  contained  energy  from  surface  forcing  which  was 
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not  considered  in  constructing  the  theoretical  decay  curves.  The  lowest  current  meters 
at  P2  and  P3  were  not  deep  enough  to  detect  the  strongest  currents  expected  near  the 
bottom.  At  P3  the  small  trapping  scales  made  it  unlikely  that  any  signal  from  a  bottom 
trapped  wave  could  be  detected.  The  current  meter  at  1000  m  was  well  above  the  e- 
folding  depth  of  diurnal  and  semidiurnal  waves.  Additionally,  the  location  of  P3  in  an 
area  with  relatively  weak  density  gradients  restricted  most  bottom  trapping  to  subtidal 
frequencies. 
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VI.     CONCLUSIONS  AND  RECOMMENDATIONS 

The  analysis  of  current  meter  records  from  moorings  located  off  Point  Sur, 
California  showed  an  increase  of  energy  at  tidal  frequencies  with  depth.  This  can  po- 
tentially be  explained  using  the  theory  of  Rhines  (1970)  for  bottom  trapped  waves  over 
the  continental  slope.  The  increased  energy  was  observed  whenever  bottom  slope,  water 
stratification  and  frequency  met  prescribed  conditions. 

Seventeen  months  of  current  meter  data  was  examined  from  mooring  P2.  The  data 
were  examined  using  complex  demodulation  and  spectral  analysis.  At  mooring  P2, 
bottom  trapped  energy  at  the  diurnal  frequency  was  observed  for  each  month  analyzed. 
Bottom  trapping  occurred  at  semidiurnal  frequencies  for  fourteen  of  the  seventeen 
months  but  the  energy  contribution  due  to  free  internal  waves,  which  influenced  the 
semidiurnal  results,  was  not  determined.  Bottom  trapping  at  mooring  P2  occurred  dur- 
ing each  month  of  the  year  and  no  seasonal  patterns  could  be  distinguished.  The  maxi- 
mum trapping  frequency  at  P2,  computed  using  N  sin  a,  was  an  accurate  indicator  of 
bottom  trapping  at  tidal  frequencies. 

No  bottom  trapping  was  observed  at  mooring  P3  using  data  collected  over  six 
months.  Weak  density  gradients  near  the  sea  floor  often  kept  the  maximum  trapped 
frequency  below  the  tidal  frequencies.  The  current  meters  at  this  mooring  were  posi- 
tioned too  far  above  the  sea  floor  for  detection  of  bottom  trapped  waves  at  diurnal  or 
semidiurnal  frequencies. 

The  bottom  trapped  wave  theory  described  by  Rhines  (1970)  was  also  used  to  cal- 
culate vertical  decay  scales  and  energy  distribution  profiles.  The  theoretical  decay  was 
modeled  using  hydrographic  data  collected  near  the  mooring  locations.    In  most  cases 
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the  modeled  decay  agreed  with  the  observed  decay  of  tidal  energy  away  from  the  sea 
floor  within  the  error  of  the  calculation. 

The  tidal  currents  off  Point  Sur  are  the  most  energetic  currents  found  over  the  con- 
tinental slope.  The  semidiurnal  currents  dominate  but  diurnal  currents  can  be  signif- 
icant, particularly  when  amplification  due  to  bottom  trapping  occurs.  Further  study  of 
bottom  trapping  should  be  conducted  using  the  continuing  current  meter  arrays  at  or 
near  the  location  of  P2.  Additional  deep  instrumentation,  particularly  just  above  the  sea 
floor,  is  needed  to  adequately  study  this  phenomena.  Because  trapping  was  observed 
throughout  the  year,  a  short  deployment  of  a  few  months  length  may  collect  sufficient 
data  to  verify  trapping  at  tidal  frequencies.  Tidal  currents  are  at  the  upper  limit  of  the 
range  of  trapped  waves  and  longer  deployments  could  be  used  to  examine  coastally 
trapped  waves  occurring  at  much  lower  frequencies  under  the  same  conditions. 
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